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Effect of vortices on the spin-flip lifetime of atoms in supercon- 
ducting atom-chips 
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PACS 37 . 10 . Gh - Atom traps and guides 

PACS 42.50.Ct - Quantum description of interaction of light and matter; related experiments 
PACS 47 . 32 . C- - Vortex dynamics 

Abstract. - We study theoretically the lifetime of magnetically trapped atoms in the close vicinity 
of a type-II superconducting surface, in the context of superconducting atom-chips. We account 
for the magnetic noise created at the cloud position by the vortices present in the superconductor 
and give specific results for our experiment which uses a niobium film. Our main result is that 
atom losses are dominated by the presence of vortices. They remain however dramatically smaller 
than in equivalent room-temperature atom-chip setups using normal metals. 



Atom chips allow to trap ultracold atomic gases in the 
vicinity of micron-sized current carrying wires [1] or per- 
manent magnetic structures [2]. Microfabrication tech- 
niques allow to design complex trapping potentials and to 
realize versatile manipulation of atoms thanks to the con- 
trol of currents or radiofrequency fields in the vicinity of 
the trapped cloud [3,4]. Atoms chips are now considered 
as a powerful toolbox that can be used for fundamental 
studies [5,6], atomic interferometry [7,8] or quantum gate 
implementation [9]. 

In many such experiments, atoms are required to be 
very close to the surface of the trapping structures. Un- 
fortunately, additional losses from the trap are experi- 
mentally observed in these conditions [10,11]. Johnson- 
Nyquist current noise in the trapping metallic wires pro- 
duce magnetic field fluctuations at the position of the 
atoms, which can induce Zeeman transitions towards un- 
trapped magnetic sublevels. This phenomenon is strongly 
enhanced in the near-field of conductors for typical spin- 
flip radiofrequencies (in the MHz range) [12,13]. The typ- 
ical geometry of these experiments is presented in Fig. 1. 
A possibility to overcome these difficulties consist in using 
cryogenic atom-chips made of superconducting materials, 
for which dissipation at RF frequencies, and hence fluctua- 
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tions, are dramatically reduced [14]. Successfull operation 
of superconducting atom-chips has been reported [15-17], 
with the aim of developping new hybrid atomic-solid-state 
systems. Concurrently, theoretical studies have made 
quantitative prediction for the lifetime increase with re- 
spect to normal metals. However they strongly depend 
on the model of superconductivity which is used in the 
calculation. Most recent articles agree to predict an en- 
hancement of at least 6 orders of magnitude [18-20]. 

In this Letter, we emphasize the importance of the 
vortex dynamics in the superconducting material on the 
atomic losses. In current atom-chip experiments, DC mag- 
netic fields of the order of 10-100 G are applied orthogo- 
nally to thin layers of type-II superconductors. Hence, we 
expect the thin film to be in the mixed-state phase, with 
vortices present in the superconducting material. Ref. [21] 
showed that the random hopping of a vortex line from a 
pinning site to another could affect the trap lifetime. We 
stress here that, in addition, the motion of the vortex line 
submitted to a RF field is in itself responsible for dissi- 
pation, a phenomenon known as "flux-flow" . However, 
this motion is partially suppressed because of the pin- 
ning of the vortices on material defects. Typical spin- flip 
frequencies are significantly smaller than vortex pinning 
characteristic energies (in the GHz range). Many theoret- 
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Fig. 1: An atom in an initial state \i) is trapped at position 
(0, 0, d) in vacuum near a film of metal or superconductor of 
thickness h. We assume that the substrate under the film is 
dielectric and does not affect the atomic lifetime. At the level 
of the atomic cloud, an external magnetic field sets the quan- 
tization axis along the x direction. 



ical and experimental studies have already been carried 
out on the dissipation of type-II superconductors in this 
low-frequency regime. We note that observations are ac- 
counted for only if one assumes a non-local response of the 
material to an applied electromagnetic field. We present 
in this letter an adaptation of the theoretical framework 
developed in Rcf. [12] to this particular non-local situa- 
tion. We have adapted the theory describing the vortex- 
dynamics in niobium [22-24] to the situation of atom- 
chips. On the basis of the measurements of the previous 
references, we evaluate quantitatively the influence of vor- 
tex dissipation on atomic lifetime for the particular case 
of our superconducting atom-chip experiment [15, 17]. 

Spin-flip lifetime calculation. Let us first con- 
sider the simple case of an infinite thickness conducting 
slab [h — > oo). The atom can decay towards an untrapped 
state [/), u> being the frequency of the i —* f transition. 
We define ko — uj/c. The contribution of a semi- infinite 
space to the spin-flip rate can be calculated in term of the 
field Green's functions [12] which is equivalent to evalu- 
ating the field radiated by the atom onto itself [25]. This 
field can be decomposed into propagating and evanescent 
plane waves (Weyl decomposition). Each of these waves 
is reflected by the surface according to Fresnel laws before 
going back to the atom. One obtains [18]: 
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where T^{w) = ^o(A*sffs) fco/(247rft) is the spin-flip rate 

in vacuum, n t h = l/(e fc s T — 1) is the mean photon number 
at frequency ui, [Lb is the Bohr magneton and gs the gy- 
romagnetic factor of the electron. The integration factor q 
is such that qko is the modulus of the wave vector compo- 
nent parallel to the surface. Evanescent waves correspond 
to q > 1. 



If we now consider the case a material described by 
a local dielectric permittivity e(u>), the polarization- 
dependent Fresnel coefficients are: 



r a (q) = 



770(g) - v(u,q) 



r P {q) = 



e(u)r) Q (q) - rj(uj,q) 
e(uj)r)o(q) +r)(w,q) 

(2) 

where rjo(q) = yl — q 2 and rj(q,uj) = y 1 's(ui) — q 2 . For 
a metal described by the Drudc model, the permittivity, 
much larger than 1, is related to the local conductivity a: 
e(u) = 1 + (equj) w ia/(eou>). The characteristic length 
associated with the material response is the skin-depth 5 = 
s/2/ (poow), typically in the (i.m range for good conductors 
at rf frequencies. The semi-infinite slab assumption holds 
then for h ^> 5. 

In the near field regime, we have l/fco ^> d. Moreover if 
we make the experimentally reasonable assumption d^$> 5, 
the main contribution to the integral of Eq. (1) is for q 
values such that l<g< y|e(w)|. One thus obtains an 
asymptotic expression of Eq. (1) with 770 ~ iq and 77(a)) w 
\/e(uj) leading to an analytical expression for Tif [14]: 
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For the case of superconductors, different theoretical 
models for the material response have been used in order 
to evaluate the spin- flip rate [18-20]: the phenomenologi- 
cal two-fluid model, the BCS microscopic model, and the 
Eliashberg theory which takes into account the scattering 
of Cooper pairs by the phonons. All those models predict 
a local complex conductivity a = o\ + ioi. As soon as the 
temperature T is significantly below the critical tempera- 
ture T c , one has cr 2 3> <J\. Equation (3) then becomes [18]: 
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The spin-flip rate is then reduced by more than 6 orders 
of magnitude as compared to the case of normal metal at 
similar temperature. 

Adaptation to the case of a type-II superconduc- 
tor. — We turn now to the situation where vortices are 
present in the film. In order to find a relation between 
the vortex dynamics and the electromagnetic radiation, 
we assume that the atom-surface distance is much larger 
than the intervortex distance ao- In this situation, the 
response of the vortex lattice to an electromagnetic field 
can be treated like that of a continuous complex hydro- 
dynamic system, taking into account vortex pinning as 
well as vortex interactions [26-28]. The theoretical mod- 
els that consider a local response of the mixed-state to 
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the electromagnetic field [29,30] underestimate the vortex 
dissipation at low frequency. It is necessary to consider 
a non-local response of the superconductor [22,23] which 
makes it impossible to define a local dielectric constant 
or a local conductivity for describing the material. Hence 
Eqs. (2-5) do not apply directly. 

In the case of a non-local response of the supercon- 
ductor, dissipation is well described in term of surface 
impedance Zs = poEs/Bs at frequency w, where Es 
and Bs are the tangential electric and magnetic fields on 
the surface respectively. The use of Zs allows to include 
the detailed microscopic response of the superconducting 
medium into a linear local relation between the tangen- 
tial electric field at the surface and the surface current 
K = Es/Zs- The rate I\/ being related to the dissipa- 
tion in the material, we expect it to be proportional to 
Re(Z s ). 

In order to find the relation between Tif and Zs, 
we express the Fresnel coefficients in terms of surface 
impedance: 



rt^Zs + Z^ 
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where Z — po c is the vacuum impedance. For the radia- 
tion of a dipole at a distance d the wave vector amplitude 
koq values mainly contributing to the integral in Eq. 1 are 
of the order of d^ 1 . For all reasonable superconductor 
models and for the range of distances considered above 
(d 3>1 |J.m), these wave vectors are much smaller than the 
inverse of the field penetration depth at the frequency u. 
The surface impedance Zs is then independent of q and 
equal to its value at normal incidence. Substituting Eq. 6 
into Eq. 1 and performing the integration is equivalent to 
the substitution e(w) = Zq/Zs in Eq. 3: 
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Note that that if we replace the surface impedance by 
its standard value in the case of a normal metal Z™ et = 
(I — i)/(aS), we exactly recover Eq. (4). 

Let us stress that, within our approximations, the sur- 
face impedance Zs allows to calculate exactly the field 
radiated by the atom in the z > half-space and hence 
the transition rate. Equation (7) is thus valid even for a 
finite thickness slab. 

Two-mode non-local response of the vortex lat- 
tice. — In order to derive Zs for a type-II supercon- 
ducting material, we evaluate the vortex response to an 
external oscillating magnetic field. As we restrict our- 
selves to the situation where the distance d is larger than 
the intervortex distance, it is possible to consider aver- 
aged macroscopic quantities for local electric and magnetic 
fields as well as for supercurrent densities. In this frame, 
the response of the superconductor can be derived from 
an equivalent of the Ginzburg-Landau free energy relat- 
ing the macroscopic quantities [26]. For this purpose, it 



is necessary to introduce the vortex field Bq = nyiptyu , 
where ~u is the local direction of the vortex lines, ny 
the density of vortices per unit area and ipo = h/2e is 
the quantum of flux. Bq is related to the macroscopic 
magnetic field B by the generalized London equation 
B + hqX^V x V s — B , where V s is the averaged macro- 
scopic velocity of the Cooper-pairs and A^ is the London 
length. In absence of vortices (B = 0), one recovers the 
London equation which leads to the Meissner effect. In 
the mixed state, a fraction B of the applied magnetic 
field penetrates the film. In the presence of a non-uniform 
supercurrent (V x V s ^ 0) vortex and magnetic field lines 
do not coincide. 

We consider now the case of an electromagnetic field 
arriving at normal incidence on the superconducting 
medium according to the geometry of Fig. (2). The vor- 
tex lattice is initially in its equilibrium position Bq = 
B = Be~l. The magnetic and vortex fields experience 
a small perturbation B(z) = Bel + b exp(ikz — iu)t)e^, 
~T?(z) = el + v a exp{ikz — iujt)e^, where \uq\ is the an- 
gle of the vortex line with the z-direction on the surface 
and \v{z)\ the same angle at the position z. Refs. [27,28] 
show that the field in the medium is a superposition of 
two propagation modes with wave vector k / and ky . The 
characteristic length of penetration Sf associated to the 
first mode is related to the flux-flow dissipation when a 
vortex line is moving: 



l 2p N B 
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B C2 is the second critical magnetic field of the material, 
and pn is the resistivity of normal electrons in the mate- 
rial at temperature T. As a consequence Sf is of the or- 
der of typical skin-depths for metals at low temperature. 
We define the complex penetration length A/ = 1/ikf = 
5f(l+i)/2. 

The characteristic penetration length for the second 
mode is: 

Av = ^ = A M /^2fLl. (9) 
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where e/ is the vortex line potential [26]. The quantity 
ifoei is the energy that is required to enter one unit of 
length of vortex into the medium. The penetration length 
Ay is of the order of Ai,. In usual superconductors one has 
Ay <§C Sf. The second mode describes the non dissipative- 
screening of the magnetic field by supercurrents located 
close to the surface. 

The incoming field must be decomposed onto these 
two modes of propagation. For each mode we can de- 
fine the displacement of the vortex line at the surface 
u j=(f,v) = J^oo v i \z)dz. The relative weight of each mode 
is determined by the boundary conditions at the surface 
for the magnetic and vortex fields. It strongly depends on 
the surface geometry of the sample. We present in Fig. 2 
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Fig. 2: Amplitude of the transverse magnetic field 6 (solid line) 
and displacement of a vortex line u (dashed line): (a) in the 
case of a flat surface, (b) in the case of a rough surface with 
characteristic length I. For both situations, the vortex direction 
v must end perpendicularly to the surface. 



two different cases of propagation in the superconduct- 
ing medium for a perfectly flat (a) and a rough surface 
(b). In both cases the vortex field must end perpendicu- 
larly to the surface. In the second case its displacement 
is prevented by its pinning on surface defects. It was pro- 
posed in Refs. [22, 27] to link the vortex displacement at 
the surface uq = Uf + uv to the angle v$ by the relation 
uq + Ivq = 0, where I is the phenomenological slippage 
length that characterizes the material. Following Fig 2(b), 
one clearly sees that I is related to the roughness of the 
surface but it also takes into account the interaction be- 
tween vortices that forces a collective response of the whole 
lattice [24]. 

The complete determination of the surface magnetic and 
electric fields, together with the vortex displacement is 
presented in Refs. [27, 28] for a infinite half-space. The 
resulting surface impedance is: 



B(l + X v )X f 
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(10) 



Finite-thickness effects. — As shown below, in our 
experimental conditions, the slab thickness h =1 |a.m is 
of the order of or smaller than Sf. The superconducting 
medium thus cannot be described as a half-space and finite 
size effects have to be taken into account. In the case of 
a type-II superconductor, they have been studied both 
theoretically and experimentally [23], in a regime where 
the magnetic field is the same on each side of the slab. It 
does not correspond to our experimental situation because 
the presence of the magnetic dipole breaks the symmetry 
between the two sides of the slab. 

We have calculated the surface impedance Z$ in the 
case of the reflection of an incident electromagnetic wave 
arriving at normal incidence on a superconducting slab of 
finite thickness. It requires to take into account 7 modes of 
propagation: the incident, reflected and transmitted fields 
as well as the evanescent propagation modes in the film 
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Fig. 3: (a) Electric (solid line) and magnetic (dashed line) field 
amplitudes in the superconducting slab of thickness h =1 u.m. 
The superconductor characteristic parameters Sf, Ay and I cor- 
respond to the niobium layer used in our experiment (see sec- 
tion "Numerical results"). The bias field applied perpendicu- 
larly to the slab is B = 100 G. The fields are expressed in units 
of the incident electromagnetic field amplitude, (b) Real part 
of the surface impedance as a function of h (in units of Sf) for 
different external bias field B 



with wave vectors ±kf and ±fey. Figure 3(a) presents the 
analytical results for the variation of the electric and mag- 
netic fields inside the superconducting film. We observe 
that only a very small part, of the order of 10 -8 , of the in- 
cident field penetrates the slab, as expected for any metal. 
Moreover, most of the magnetic field is screened by super- 
conducting currents carried by the the modes ±fcy which 
do not dissipate. Hence an even smaller fraction of the in- 
coming wave is dissipated by the vortex displacement and 
contributes to Re(Zs). 

We present in Fig. 3(b) the calculated real part of Zs 
as a function of the slab thickness for different external 
field values B. We recover the value Re(Z^) of Eq. (10) 
as soon as h > 25 f. In our experiment, we are in the 
opposite limit h <C Sf, for which we fit the results with 
the phenomenological formula: 



Re[Z s (h)] « ~Re(Z?) 



(11) 



Numerical results. — In the specific case of our su- 
perconducting atom-chip [15, 17], we can evaluate all rel- 
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evant parameters and determine Zg- Resistance mea- 
surements of the sputtered Nb film give a normal re- 
sistivity above transition to the superconducting state 
Pat =15 cm. The comparison with measurements at 
room temperature gives a residual resistance ratio (RRR) 
of 4.6, which indicates that the film is in the so-called 
"dirty limit" . Measurements of the critical magnetic field 
for similar films in this limit gives B C2 — 4.5 • 10 4 G, a fac- 
tor 15 larger than the pure case value [31]. To first order 
the product B Cl B C2 is a quantity that weakly depends on 
the quality of the film and remains almost constant [32]. 
Hence, we expect B Cl to be a factor 15 smaller than the 
pure case value [33], i.e. B Cl = 80 G. The vortex po- 
tential e; depends on the applied magnetic field B. For 
normal operation of an atom-chip, we fulfill the condition 
B ;=y B Cl <C B C2 and we have fi a ei w 0.9B Cl w 70 G [34]. 

A microscopic calculation of I is out of the scope of this 
Letter. It can be found in [24]. Nevertheless, it is possi- 
ble to evaluate this quantity by linking it to the critical 
current in the superconducting slab I c — JJ j s d r. On 
the one hand, we have I c sa 2weiv c [23], where w 3> h 
is the width of the slab and v c is the maximum angle at 
which the vortex line can bend before their pinning on the 
surface breaks. On the other hand, the critical angle is 
reached when the displacement u is of the order of the in- 
tervortex separation an = \J ifio/B. Combining those two 
equations with the boundary condition between u and v , 
we obtain: 




We will now consider that the atoms are trapped above 
the ft. =1 urn thick superconducting Z-wirc of Rcf. [15]. 
We consider this finite width conductor (w =40 |0.m) as 
an infinite film, leading to a worst case estimate for the 
spin-flip rate. We assume that an homogeneous external 
bias field B = 100 G is applied perpendicularly to the 
slab. This crude assumption corresponds to a worst case 
situation. Another external magnetic field parallel to the 
surface determines the spin-flip frequency uo = 2tt x 2 MHz. 
Eq. (8) gives then Sf =9.2 urn and using Eq. (9) we obtain 
Ay = 29 nm (with A^ = 45 nm [33]). In presence of a cur- 
rent in the slab, the value I can be obtained by replacing 
Ic by I c — I in Eq. (12) where I = 1.4 A is the actual 
current in the slab required for trapping (7 C = 1.76 A). 
We get I — 250 nm. Figure 4 compares the spin-flip 
lifetime t = T^y derived from Eqs. (11) and (7) as a 
function of the distance d. These predictions are com- 
pared to our measurements for a gold layer of thickness 
h = 200 nm [35]. We also present the results of the BCS 
model of Refs. [18-20]. 

Our predictions stand in a regime where the distance 
d 3> Sf, Ay, an. The results of Fig. 4 are therefore valid 
for d >10 um. In this regime the vortex dissipation is the 
limiting factor for the spin-flip lifetime. We find a reduc- 
tion of about 3 orders of magnitude of the lifetime with 
respect to the two-fluid model. However, superconducting 
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Fig. 4: Spin-flip lifetime as a function of atom-surface dis- 
tance as predicted by Eq. (7) (solid line). It takes into ac- 
count the measurements of our real superconducting film prop- 
erties and its finite thickness h —1 [im. The dashed lines cor- 
respond to two ideal cases where the thickness of the film is in- 
finite (long dashes) , or its superconducting properties are ideal 
(clean limit, short dashes). The dotted line corresponds to 
the superconducting BCS model [19]. The experimental points 
correspond to our measurements for a Gold layer of thickness 
h — 200 nm [35]. They are in good agreement with the pre- 
dictions of Eq. (4), with (TAu,4.2K = 6.7 x 10 9 fi m(dot-dashed 
line) 

Nb remains significantly better than normal metals. 

We have also represented in Fig. 4 the predicted life- 
times in the case of an semi-infinite superconducting film, 
where the surface impedance is given by Eq. (10), or for 
a film with ideal purity and surface quality (clean limit). 
It is interesting to note that these two situations corre- 
spond to a degraded lifetime compared to the dirty film. 
This might seem paradoxical. In this clean superconductor 
regime however, surface pinning is significantly reduced, 
the vortices can move with a larger amplitude. In addi- 
tion the change of the normal fluid conductivity increases 
the viscous drag of the vortex lattice. Two phenomcnons 
therefore contribute to the increase of the dissipation. 

We also note that we have considered here a worst case 
limit where we assume that all the surface is subjected 
to a an external field of 100 G and a homogeneously dis- 
tributed current 1.4 A. A more precise calculation of the 
the atomic losses should include the inhomogeneous cur- 
rent and vortices distribution on the surface. Neverthe- 
less, our results show that vortices play a crucial role in 
the spin-flip lifetime of atoms in the close vicinity of a 
type-II atom-chip. 

In conclusion we have adapted the formalism of atomic 
spin-flip lifetime of an atom close to a metallic surface to 
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the non-local elcctrodynamic response of the vortex lattice 
in type-II superconductors. Note that this model should 
also described type-I superconductors, as a thin film of 
such a material will contain vortices. On the other hand, 
lifetime close to superconducting materials remains signif- 
icantly better than close to normal metals. Our results 
predict a lifetime of 10000 s at 20 urn from the surface 
opening new perspectives for the coherent manipulation 
of ultracold atoms in the vicinity of superconductors. 
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